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ABSTRACT 

We report the detection of the first 2 microlensing candidates from the Wendelstein Calar Alto Pixellensing 
Project (WeCAPP). Both are detected with a high signal-to-noise-ratio and were filtered out from 4.5 mill, pixel 
light curves using a variety of selection criteria. Here we only consider well-sampled events with timescales 
of Id < ffwhm < 20 d, high amplitude, and low of the microlensing fit. The two-color photometry (/?,/) 
shows that the events are achromatic and that giant stars with colors of [R-I] ~ 1.1 mag in the bulge of M31 
have been lensed. The magnification factors are 64 and 10 which are obtained for typical giant luminosities 
of M/ = -2.5 mag. Both lensing events lasted for only a few days (ff^^m ~ ^-^^^ ^'^'^ ^^hm ~ The event 

GLl is likely identical with PA-00-S3 reported by the POINT-AGAPE project. Our calculations favor in both 
cases the possibility that MACHOs in the halo of M31 caused the lensing events. The most probable masses, 
0.08 Mq for GLl and 0.02M3 for GL2, are in the range of the brown dwarf limit of hydrogen burning. Solar 
mass objects are a factor of two less likely. 

Subject headings: gravitational lensing — galaxies individual: M31 — dark matter — galaxies: halos 



1. INTRODUCTION 

Microlensing experiments are an ideal method to search for 
dark objects within and between galaxies. A large number of 
microlensing events have been detected towards the Galactic 
bulge constraining the number density of faint stars in this di- 
rectionjAfonso et al. 1999; Alard 1999: Alcock et al. 2000a; 
lUdalski et a l. 2000). Towards the LMC only 13-17 microlens- 
ing events have been reported so far (lAlcock eral.n20001? ). 
If all this events are attributed to 0.5 M© MACHOs, the as- 
sociated population of dark objects would contribute up to 
the 20% level to the dark matter content of the Milky Way 
(|Alcock et al. 2000b). However, both the relatively large size 
of the LMC relative to its distance and the nature of the lenses 
has cast doubt on this interpretation. It is indeed likely that a 
large fraction of the microlensing events towards the LMC are 
due to self-lensing of stars within the LMC ( Evans & Kerins 
(|2000); Lasserre et al. (2000) and references therein). 

Studying microlensing events towards M31 allows to sep- 
arate self-lensing and halo-lensing in a statistical way, since 
the optical depth for halo lensing is larger on the far side of 
M31. In M31 individual stars can not be resolved and one 
therefore has to us e the pixellensing technique ( Baillon et al. 
ll993^ICrottsll992^ to follow the vai-iability of sources blended 
with thousands of other sources within the same pixel. First 
detections of possible microlensing events were reported 
by several pixellensing experiments ( Ansari et al. _199a 
lAuriereetal. 2001; Calchi Novati et al. 2003t ICrotts et all 
[T99fc Paulin-Henriksson et aL .200Z .2003.) . But since the 
candidate nature of only 5 of these events is convincing, no 
conclusions concerning the near-far asymmetry or the most 
likely dark matter lensing masses could be drawn yet. 

The Wendelstein Calar Alto Pixellensing Project (WeCAPP, 
iRiffeser etaDl200l started in 1997 with test observations. 

' Based on observations at the Wendelstein Observatoi'y of the Univer- 
sity Observatory Munich and at the German-Spanish Astronomical Center, 
Calar Alto, operated by the Max-Planck-Institut fiir Astronomic, Heidelberg, 
jointly with the Spanish National Commission for Astronomy. 

^ arri@usm.uni-muenchen.de 



Since 1999 the bulge of M31 was monitored continuously 
during the time of visibility of M3 1 . The analysis of our 4 year 
data will allow not only the identification of very short dura- 
tion events (eg., in the 4* year data of the combined field have 
been taken on 83 % of possible nights) but also the separation 
of long duration (< 20 d) microlensing events from long pe- 
riodic variables as Mira stars. For this letter we analyzed the 
short duration events (ffwhm < 20 d) within one season of Calar 
Alto data and restricted the detection to high-signal-to noise, 
high-magnification events. We report our first 2 microlensing 
candidates of that type. 

2. OBSERVATIONS AND DATA REDUCTION 

WeCAPP monitors the central region of M31 in a 17.2 x 
17.2 arcmin^ field with the 1 .23 m telescope of the Calar Alto 
Observatory. In addition, a quarter of this field, pointing to- 
wards the far side of the M3 1 disk along the SE minor axis, 
was observed with the 0.8 m telescope of the Wendelstein Ob- 
servatory. The data analysis and candidate selection reported 
in this letter is based on the season from June, 23th 2000 to 
February, 25th 2001 and is restricted to the Calar Alto data 
only. During this period, M31 was observed during 43% per- 
cent of all nights. Observations were carried out in R and / 
filters close to the Kron-Cousins system. We estimate the sys- 
tematic error in the (R-I) color to be < 0.05 mag. 

We have developed a pipeline based on Gossl & Riffesed 
( I2OO2I) and lRiffeser et al] ( fin prep. I) , which performs the stan- 
dard CCD reduction, position alignment, photometric align- 
ment, stacki ng of frames , PSF matching using Optical Image 
Subtraction (lAlardll998h . and the generation of difference im- 
ages. For the data presented here all data within one night are 
coadded, yielding one difference image per night. The re- 
duction package includes full error propagation for each pixel 
through all reduction steps. In this way, all data points are 
properly taken into account in the search for variables. 

3. SELECTION CRITERIA 

We investigate only pixels which have more than 10 data 
points in R and /, which applies for 85% of the 2K x 2K 



2 



Riffeser et al. 



Table 1. Selection criteria 



selection criterion 



number 



analyzed liglit curves 4492250 

light curves with > 10 data points 3835407 

3 successive 3-cr in R or / 5 17052 

x| < 1.3 andx? < 1.3 186039 

1 day < rfwhm < 20 days 9497 

3-(T light point inside ffwhm 1829 

sampling: sidei > 20%, side2 > 5% 256 

feff > 10 mediaUenw in R and / 15 

candidates 2 



field. For each pixel we define a flux baseline by iterative 3a- 
clipping of all outliers with higher flux. All pixels which have 
at least 3 successive (positive) S-cr deviations from this base- 
line are considered as variables. We fit the m icrolensing light 
curve for high-amplification events (>G ouldlll996.) simultane- 
ously to the R- and /-band pixel light curves for every vari- 
able. The fit has 6 free parameters: full-width-half-maximum 
ffwhm and the time to of maximum amplification (these 2 pa- 
rameters are the same for both filters), amplitude Feff,« , color 
Feffj / Ff-ffK , and baseline levels cr and cj. Variables with a 
reduced x| > 1.3 or xj > 1-3 are discarded. In this way we 
exclude light curves that are not achromatic or not symmet- 
ric. We also exclude events with ffwhm > 20 days which can 
be confused with long periodic variables like Mira stars, as 
long as only one season of data is investigated. In addition, 
all candidates which do not have at least one significant data 
point (3 (7 deviation from the baseline) within ffwhm of the time 
of maximum amplification are rejected. We further define the 
sampling quality for the falling and rising parts of each light 
curve within {Iq- 15d,fo) and (fo,fo + 15d): within these time 
intervals we require a sampling of the area under the light 
curve of at least 20% on one side and of at least 5% on the 
other side (Table0. 

Here, we present only the two microlensing candidates 
which have amplitudes 10 times larger than the median error 
of the light curve (see Figure[0. Both candidates fit perfectly 
to a symmetric microlensing light curve. Ruling out system- 
atic offsets for the points and errors on the trailing side of GL2 
(which is strongly proved by the 6 single images of that night 
in each filter) a non microlensing light curve of a variable 
source hardly fits the data points of GL2. Both microlens- 
ing candidates are detected in several pixels (11 for GLl and 
4 f^or GL2) inside the PSF of the position of the lensed object. 
This explains the reduction from 15 events to 2 events in the 
last line of Table [0 The amplification light curves were ob- 
tained by calculating the total flux within the PSF area of each 
microlensing event. For both candidates, the selection criteria 
exclude variable stars like Miras, novae or dwarf novae. Ex- 
tracting lensing events with less good time sampling, lower 
amplitude, or events located close to other variables requires 
refined selection criteria. These will be discussed in a future 
paper which will also include a test of the detection efficiency 
and false detection rate with Monte Carlo simulations. 

4. MICROLENSING CANDIDATES 

The parameters of both lensing candidates are summarized 
in Table 121 Their light curves are shown in Figure GLl, 
the highest S/N lensing event candidate in our sample, lies 
4"1 to the SW of the nucleus of M31. GL2 is 4"4 to 
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Fig. 1.— Light curves of WeCAPP-GLl and WeCAPP-GL2. The I-band 
light curve (red symbols, right axis) has been scaled to the R-band light curve 
(blue symbols, left axis). The scaling factors were deiived from the lensing 
fit (black curve) and con'espond to a color {R—I) of 1.05 for GLl and 1.08 
for GL2. In addition we show the r' and i' data from the POINT- AG APE 
PA-00-S3 event (green symbols) scaled to our data. 



Table 2. Positions and parameters of the microlensing candidates 



Name 


WeCAPP-GLl 


GLl & PA-00-S3 ' 


WeCAPP-GL2 


02000 


00*42"' 30:^3 


00*42'"30?3 


00''42"'32?8 


(52000 


+41°13'00"8 


+41°13'00:'8 


+41°19'56"5 


to [JD] 


2451850.80±0.13 


245 1850.84 ±0.02 


2451894.77±0.2I 


'fwhm [d] 


1.38±0.53 


1.65 ±0.10 


5.41 ±2.49 




13.4 ±5.4 


12.4±0.6 


1.7 ±0.5 




28.0± 11.2 


25.7± 1.5 


3.6±1.1 


(R-I) 


1.05 ±0.08 ±0.05 


1.05 ±0.08 ±0.05 


1.08 ±0.24 ±0.05 




1.23 


1.22 


1.02 



'deiived from a fit to the total set of data points (WeCAPP & POINT- 
AGAPE) 



the NW of the nucleus. Our data have been astrometrically 
cal ibrated using bright f oreground stars observed with HST 
by Jablonkaet al. ( 1999) and ground based observations by 
Magnier et al. ( 1992). Our two calibrations agree within 0.5" 
in declination and 0.7" in right ascension, well consistent 
with the astrometric accuracy of 0.8" to 1.0" of the Mag- 
nier et al. catalog. After we had detected GLl and GL2, 
we cross-checked with events reported by the POINT- AGAPE 
survey for the same perio d of time and the same field in M3 1 
( Pauli n-Henriksson et al.lf2003iD . It appears that GLl is likely 
identical with PA-00-S3 which occurred at the same time 
(Figure [Hi. Because POINT- AGAPE did not provide a flux 
calibration of their data, we had to assume a scaling factor for 
the amplitude. The zero point in time was not adjusted. The 
data points from WeCAPP and POINT- AGAPE complement 
each other nicely and make GLl the best pixellensing event 
found so far in M3 1 . GL2 also falls in the observing period 
covered by POINT-AGAPE but their time sampling around 
the event is poor. This may be the reason why GL2 was not 
detected. 

The parameters of t he lens ing fit are degenerate for high 
magnifications (Gouldl [l996h . i.e. for amplitudes Aq ^ 1 
which correspond to impact angles much smaller than the Ein- 
stein angle 9e. In this case, we obtain for the light curve: 



Fo[A{t)-l]~^F,s 
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^fwhm' 



■(f-f0)2 + l 
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Table 3. Amplifications for different source star luminosities 



iVll 


Ao 




GLl 






-1.9 


120 


47.4±18.1 


-2.5 


69 


27.2 ±10.4 


-3.7 


24 


9.0±3.5 


1.1 


8.2-10^"^ 


(3.3±1.3)-105 


GLl 


& PA-00-S3 




-1.9 


110 


52.1±3.2 


-2.5 


64 


30.0±1.8 


-3.7 


22 


9.9 ±0.6 


1.1 


7.6 - 10^ 


(3.6±0.2)-105 


GL2 




-1.9 


16 


23.9 ±11.0 


-2.5 


10 


13.8±6.3 


-3.7 


4 


4.6±2.1 


8.3 


1.8 - 10^ 


(2.9±1.3)-105 



with intrinsic flux of the source Fq and F^^ = Fq/ uq w FqAq, 
where uq is the impact angle in units of the Einstein angle. 

In order to obtain Aq we evidently need to know the source 
flux Fq. We can get rough constraints by considering the 
color of the light curves which due to our selection criteria 
are achromatic. We obtain {R-I)gli = 1-05 and {R-I)gl2 — 
1.08. For main sequence stars this converts^ iCassisi et alJ 
Il998h into absolute magnitudes of M/ w 8 mag. If the sources 
are giants, then the magnitudes are Mj ~ -2.5 mag with 
a minimum of -1.9 and maximum of -3.7 ( lBesselllll979t 
lGri11mairetaljT996i). 

We now derive the lensing parameters as follows: 

Fq,I ffwhm (vtfE)^c2Ds 
"0 = 'E^ > — 7= > Miens — FTT 

where is the Einstein timescale, v, is the transverse veloc- 
ity between source and lens, and Di and are the angular 
distances to the lens and source. 

To estimate plausible lens masses the Einstein timescales fe 
are calculated for fixed luminosities of possible source stars 
(see Table|3}. Note that the errors in reflect the accuracy of 
the determination of ffwhm in the degenerate Gould fit only, 
and do not account for the systematic uncertainties due to 
the unknown luminosities of the sources. If the source is a 
main sequence star, we need very high magnifications: typi- 
cally Ao = 8 • 10^ and Ao = 2 • 10^ for GLl and GL2, respec- 
tively. The corresponding lens masses (for v, ~ 210km/s, 
Dl = 768 kpc, Ds = 770 kpc) are M w IO^Mq, an implausi- 
bly large value. If the source is a giant, the required magnifi- 
cations are reduced to Aq.gli = 64 and Ao.gl2 = 10 (see Table 
13; the typical self-lensing masses become M ~ 0.8 Mq and 
M — O.2M0, i.e. they are typical for low mass stars. 

Assuming the source to be a red giant with M/ = -2.5 mag 
we calculate the probability p{M,tE) that a microlensing 
event of observed timescale f£ can be produ ced by a lens of 
the mass M. Following the calculations of Lletzer & Massd 
((1994) eq. (8) and Jetzer (1994) eq. (11) we get 

p{M,tE) ^ ^M) I ps{Ds) I Pl(Ol)/ (^^) ^dDi^dDs 

with the mass function (MF) ^(M), the sources density 
Ps{Ds), the lenses density the velocity distribution 

/(vt) and the Einstein radius Re{Di^,M ,Ds)- 

' transformed on the observational plane bv lMarastonI 119981) 




0.01 0.1 1 0.01 0.1 1 
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Fig. 2. — Mass probability for GLl (left panel) and GL2 (right panel) 
for lens-source configurations: halo-bulge (filled circles), halo-disk (open 
circles), bulge-bulge (red), bulge-disk (green), disk-bulge (blue), disk-disk 
(magenta). The maximum of each curve is scaled to reflect the total prob- 
ability of a respective lens-source event relative to the case of a halo-bulge 
lensing event with the most probable MACHO mass. For example in case 
of GLl the probability for bulge-bulge lensing relative to halo-bulge lensing 
with 0.08M. lenses becomes 0.6 (maximum of red curve). A halo consist- 
ing of 0.014M;., MACHOs would have the same probability as bulge-bulge 
lensing. Note that the shapes of the distributions for bulge and disk lenses are 
strongly affected by the used mass function § {M) . 



The distribution of matter in the central part of M3 1 is based 
on the bulge model of Kent ( 1989). The disk is modeled with 
a radial scale length of 6.4 kpc and an exponential shape, and 
with a vertical scale length of 0.3 kpc and a sech^-shape. The 
halo is modeled as an isothermal sphere with a core radius of 
Tc = 2 kpc. The velocity distribution was calculated from a 
Maxwellian halo bulge and disk velocity distribution with an 
additional rotation for bulge and disk ( Kerins et al. 2001 ). 

For the bulge lenses we take the IMF as derived for the 
galactic bulge ^ M~^-^^ (ZoccaH et al. 2000). For the disk 
population we adopt a Gould I MF ^ ~ M~^-^^ with a flattening 
^ ^ M""^*' below 0.59M© ( G ould, Bahca ir& Flvnn 1997). 
Both are cut at the lower end at the hydrogen burning limit of 
O.O8M0. At the upper end the bulge MF is cut at the main se- 
quence turn-off 0.95 M0 (Maraston, priv. com.) and the disk 
MF at lOM©. The IMF for the potential MACHO population 
residing in the halo of M3 1 is of course unknown. We there- 
fore calculate the probability distribution for halos consisting 
of one mass only, i.e. taking (5-function IMF's centered on the 
lens mass ^ = i5(M- Miens )/A/iens- Moreover we assume that 
the whole dark halo of M31 consists of MACHOs. Lensing 
by Galactic halo objects has an order of magnitude smaller 
optical depth and is therefore neglected in our considerations. 

The results are shown in Figure|2l For M3 1 halo lenses the 
most probable masses are O.O8M0 for GLl and O.O2M0 for 
GL2. In the case of self-lensing the most probable masses are 
about a factor of 4 bigger Taking the most likely halo lens 
masses, the ratio of the probabilities that the lenses are part 
of the dark halo or the stellar content /:'haio/(/'buige + Pdisk) is 
1.6 for GLl and 3.3 for GL2. We conclude therefore that it 
is likely that lenses residing in the halo of M31 caused the 
events in both cases. 

5. discussion AND OUTLOOK 

We presented the first two high (S/N), short timescale 
microlensing events from WeCAPP. GLl is likely identical 
to PA-00-S3 found by POINT-AGAPE. Combining the data 
from AGAPE with ours shows that the error bar of the derived 
Einstein time scale becomes smaller by a factor of 5 compared 
to the individual error bar This demonstrates the importance 
of a good time sampling of the events. We derived the colors 
of the lensed stars, the amplification factors and likely lens 
masses for both bulge/disk self-lensing and MACHO lensing. 
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We showed that red giants are the Ukely source objects, while 
main sequence stars are highly implausible. 

Self-lensing in the bulge can only be separated from MA- 
CHO lensing statistically. Halo-lensing events show a spa- 
tial asymmetry because the optical depth for lensing events 
is higher for stars on the far side o f M31 than on its near 
side l^rotts 1992; F liri et alJin prepJ) . In contrast, bulge self- 
lensing is symmetric. 

The bulge self-lensing hypothesis yields lensing stars at 
or below the the main sequence turn-off of the M31 bulge. 
On the other hand, if the lensing events are caused by MA- 
CHOs, their masses are typically very low, most probable be- 
low O.1M0. Masses in the range of O.5M0 to IMp are more 
unlikely. 

So far, we have analyzed one observing season and re- 
stricted the lensing search to short-time, high-amplification 
events in order to avoid confusion with variable stars. The 
whole WeCAPP dataset will allow us to identify all variables 
and thus will enable a search for lower amplitude and longer 
duration microlensing events. 

Decreasing the amplitude threshold will increase the de- 
tected rate of events in two ways. As the event rate is propor- 



tional to the inverse of the minimum required magnification 
^o.min in the pixellensing regime we expect to detect more 
lensed giants. On the other hand lowering the amplification 
threshold could make it po ssible to detect al so highly ampli- 
fied main sequence stars (Han & Gould 1996) which exceed 
the evolved stars in the bulge of M3 1 by a factor of more than 
a hundred. How many more lensing events will be detected 
depends on the mass function of the lenses but we can expect 
at least a factor of a few (Fliri et al. in prep.). 

Finally, the effects of time sampling and noise properties 
of our sample on the detectability of lensing events have to 
be taken into account. Results of the modeling of these ef- 
fects for events of different durations and amplitudes using 
Monte-Carlo simulations will be presented in a future publi- 
cation. With the full dataset we expect therefore to increase 
the number of lensing events to detect the predicted asymme- 
try of MACHO lensing or to rule out a significant MACHO 
population in the halo of M3 1 . 



This work was supported by the SFB 375 Astro-Particle- 
Physics of the Deutsche Forschungsgemeinschaft. 
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